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PREFACE 


This report presents the results of studies conducted during the 
period July 19, 1969 — July 19, 1970, under NASA Research Contract 
NAS 8-21432, "Lunar Surface Engineering Properties Experiment Definition." 
This study was sponsored by the Lunar Exploration Office, NASA Head- 
quarters, and was under the technical cognizance of Dr. N. C. Costes, 

Space Science Laboratory, George C. Marshall Space Plight Center. 

The report reflects the combined effort of five faculty investiga- 
tors, a research engineer, a project manager, and eight graduate research 
assistants, representing several engineering and scientific disciplines 
pertinent to the study of lunar surface material properties. James K. 
Mitchell, Professor of Civil Engineering, served as Principal Investigator 
and was responsible for those phases of the work concerned with problems 
relating to the engineering properties of lunar soils and lunar soil 
mechanics. Co-investigators were William N. Houston, Assistant Professor 
of Civil Engineering, who was concerned with problems relating to the 
engineering properties of lunar soils; Richard E. Goodman, Associate 
Professor of Geological Engineering, who was concerned with the engineer- 
ing geology and rock mechanics aspects of the lunar surface; Paul A. 
Witherspoon, Professor of Geological Engineering, who was concerned with 
fluid conductivity of lunar surface materials in general; Franklin C. 
Hurlbut, Professor of Aeronautical Science, who was concerned with 
experimental studies on fluid conductivity of lunar surface materials; 
and D. Roger Willis, Associate Professor of Aeronautical Science, who 
conducted theoretical studies on fluid conductivity of lunar surface 
materials. Dr. Karel Drozd, Assistant Research Engineer, performed 
laboratory tests and analyses pertinent to the development of a borehole 
jack for determination of the in situ characteristics of lunar soils 
and rocks; he also helped in the design of the borehole jack. H. Turan 
Durgunoglu, H. John Hovland, Laith I. Namiq, Parabaronen Raghuraman, 

James B. Thompson', Donald D. Treadwell, C. Robert Jih, Suphon Chirapuntu, 
and Tran K. Van served as Graduate Research Assistants and carried 
out many of the studies leading to the results presented in this 
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report. Ted S. Vinson, Research Engineer, served as project inanager 
until May 1970, and contributed to studies concerned with lunar soil 
stabilization. H. John Hovland served as project manager after May 
1970, and contributed to studies concerned with soil property evaluation 
from lunar boulder tracks. 

Ultimate objectives of this project were: 

1) Assessment of lunar soil and rock property data using information 
obtained from Lunar Orbiter, Surveyor, and Apollo missions. 

2) Recommendation of both simple and sophisticated in situ testing 
techniques that would allow determination of engineering 
properties of lunar surface materials. 

3) Determination of the influence of variations in lunar surface 
conditions on the performance parameters of a lunar roving 
vehicle. 

4) Development of simple means for determining the fluid 
conductivity properties of lunar surface materials. 

5) Development of stabilization techniques for use in loose, 
unconsolidated lunar surface materials to improve the 
performance of such materials in lunar engineering application. 

The scope of specific studies conducted in satisfaction of these objectives 
is indicated by the following list of contents from the Detailed Final 
Report which is presented in four volumes. The names of the investigators 
associated with each phase of the work are indicated. 
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VOLUME IV 


Studies on Fluid Conductivity of Lunar Surface Materials 

Chapter 1 . THEORETICAL STUDIES 
P. Raghuraman and D. R. Willis 


INTRODUCTION 

Theoretical studies on the problem of developing a probe capable of 
measuring the fluid conductivity of porous media under lunar conditions 
have been pursued in two principal directions. In the first place ^ we 
have examined some aspects of the fundamental question as to whether or 
not the probe can operate in quasi steady state conditions or whether it 
will be operated while various transient phenomena are still affecting 
the flow. A reasonable estimate of the time to establish steady state 
conditions appears necessary in view of the limitations on the maximum 
in situ testing time imposed — both by gas storage limitations and by the 
time an astronaut can devote to this experiment, with this in mind, we 
have considered the effect of dead end pores, which initially contain 
essentially no gas and act for some time as sinks for mass flow, when gas 
first flows down the main (through) pores. Clearly a truly steady flow 
cannot be established before there is pressure equilibrium between the 
dead end pore and the open pore. In the next section, the influence of 
the length-to-radius ratio on the time constant for filling up the dead 
end pores in the free molecular limit is considered. 

In the second principal line of study, we have focused our efforts on 
evaluating the degree of sophistication necessary to analyse the flow 
through a porous media with vacumn as one boundary condition. We recognise 
that the flow, in general > would contain an initial region of continuum - 
flow which would undergo a gradual transition to free molecular flow. In 
an effort to obtain rough estimates, a sudden freeze model is proposed. 

The third section contains the application of this model to various other 
models of the porous media. 
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FREE MOLECULE TIME CONSTANT OF DEAD-END PORES 
Formulation and Assumptions 

The unsteady Boltzmann equation is used to find the time constant 
for filling up a dead-end pore at zero pressure. The dead-end pore is 
visualized as a straight cylinder (Fig. l-l) , of radius "a" and length "Z", 
closed at one end. The problem is posed as the calculation of the time 
constant of such a cylinder at zero pressure, separated from an infinite 
reservoir of gas by a diaphragm, which is suddenly removed at some instant 
t = 0. An approximate scheme is formulated to determine the time 
constant for filling the cylinder, and the influence of the length of 
th& cylinder on the time constant. 

The basic assumptions made in the formulation are: 

1. The flow is sufficiently rarefied that the effect of intermolecular 
collisions within the cylinder can be ignored compared to the 
effect of collisions with the boundary, i.e. , the flow is free 
molecular. 

2. The molecules that strike the walls undergo perfect accomodation. 

3. There is no accumulation or ablation at the wall. 

4. The boundary distribution functions are Maxwellian; that is, the 
infinite pocket of gas outside the cylinder has a Maxwellian 
distribution, as do the molecules reflecting from the cylinder 
wall. 

Taking the mass balance on a unit area of the wall at axial distance z 
and time,t, the number flux C^(z,t), from the cylindrical wall, is 

? (z,t) = C (z,t) + C (z,t) + C, (z,t) (4-1) 

c ec cc be 

while the number flux, ?j^(r,t), per unit area from the back end at radius, 
r, and time, t, is 
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respectively. Also, the total outflux, e (t) , from the cylinder at 
time t is 


e (t) = e (t) + e (t) , 
c 


(4-3) 


where e^ is the number flux to the cylinder exit from i. 

Each of the terms on the righthand side of Equations (4-1), (4-2), 
and (4-3) are evaluated using the ray tracing technique. The calculations 
detailed in the Appendix yield the following equations for (4-1) , (4-2) , 
and (4-3), respectively: 






(£ - z) ^ + a^ + r^ + 2ar cos0. 
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where and are the number density and temperature outside the 
cylinder, Tw is the temperature of the inside surfaces, and c is 
molecular speed scaled by 

We note that ? (z,t), (r,t) and e(t) are equal to zero for t S 0 

c b 

while in the limit of larger values of t, they all tend to a constant 
value equal to |^n^ Equations (4-4) and (4-5) constitute an 

integral equation for C^(z,t). Once it is solved, ^^^(r,t) and, hence e(t), 
can be easily determined from (4-5) and (4-6) , respectively. The time 
constant is then determined as that value of time t at which e(t) reaches 
(1 - 1/e) of its final value j^n^ ^2RT^^^J/2 tt^. 

We will concern ourselves only with those situations wherein 3.nd, 

T do not differ appreciably. Then the only characteristic time in the 
filling of the cylinder is a/^2RT^^ |, and the time constant will be 
equal to A^il/a^ ja/ ^2RT^^^J . The problem as posed is to determine the 
function A(Jl/a), which gives the dependence on (fi./a) . In view of the 
motivation of this work, it does not seem appropriate to solve Equations 
(4-3) , (4-4) , and (4-5) in detail. Rather, a series of assumptions are 
made to yield simpler equations; these are then used to get the upper and 
lower bounds for the time constant as a function of the length-to-radius 
ratio, Jl/a. In the next two subsections, the assumptions and the final 
equations for the upper and lower bounds are detailed. 

Upper Limit to the Time Constant 

The various assumptions made are: 

1. The molecular velocity C(2RT^j involved in the time of flight 
occurring in the various number fluxes, is replaced by the most 
probable velocity ^2RT^^^. 

s 2. The center of the back end is taken as a typical and representa- 
tive point for the back end. 

These two assumptions lead to simplifications of the various terms 
involved on the righthand side of Equations (4-1) , (4-2) , and (4-3) . The 
details are elaborated in the Appendix. The simplification for each term 
follows the derivation of the respective term. 


I 
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3. The lengths of flights involved in the fluxes are square roots 
of expressions involving Si, a, r, 9, and z in various combina- 
. tions. It is proposed to get rid of the square root and represent 
these' lengths by "average" values. 

Thus, in Equation (4-4), we represent j^(z - z')^ + 2a^(l + cos6)j 

averaged over 0 by |z' - z| + 2a/k, with k, a function of (| z - z'| )/a, 

and |(i!, - z) ^ + a^ + r^ + 2ar cosSj at r = 0 (due to assumption 2) by 

(£, - z) + a/k , where k is a function of (Si - z')/a. In Equation (4-5) 

2 2 1 

we represent f(il - z')^ +a^ +r^ + 2ar cos9| at r = 0 (due to assumption 2) 


[(il - z’ 


by (S, - z') + a/k where r is a function of (Si - z')/a. Finally in 
^ r 9 2 

Equation (4-6) we represent j z“ + a + 2ar cos9 + r i averaged over 6 and r 

L j- J - 

by z + 2a/kg where k^ is a function of z/a, and + r^ + r'^ + 2rr' cos0 j 

at r = 0 (due to assvimption 2) , averaged over r* , by il + a/k^ with k^^ as a 

function of Z/a. The last averaging is, however, valid only for large or 
"medium" length- to-radius ratios. For short length to radius ratios it is 
felt to be appropriate to represent the average of (X- ^ + r ' by S- itself . 


Finally, we naively assume that k^ = k '2 = k^=k|^ = a constant = k, say. 
This carries the tacit assumption that the exact values of k^ , k^, k^, and 
k^ do not have a profound bearing on the final result. 

The following scaling is introduced 
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?c(2.t) 


[ne(2RTe)V2(TT)^] = f(|j^ I, Tj) 


Cjj(0,t) 




■ t' 


> (4-7b) 




L"^ J 


The shorthand notation F(I,J), E(J) , and F (LL + 1, J) are used for 

F |cll + 1) ^ , TjJ , respectively. The scaled flux 
from the back end has been called F(LL + i, J) for easier understanding of 
the mechanics of the problem. 

If J,K, and I are taken as continous variables, then consequent to the 
three assumptions. Equations (4-4) through (4-6) take the following form; 


F(I,J) = 

TT 


^ - 1 2Nk 




d0 sin^e’’ ^ ^ 


( 


X - 


tan^e 

4N^k2 


f ^ f a, [1 l’^ - - ’^y * 

^ J L [(i - k)^ * 4„v]^ . 


X F(K, J - |k - l| - 2N) + 


(LL - I) 
2Nk 


(LL 

- I)^ + 2N^k^ 


V - ' 

(LL - I) 

^(LL - I)^ + 4N^k^] 


F ( U., I 1 , .1 - LL + 1 - N ) 


(4-H) 



with 


F(LL 


+ 1, J) 


. .-(*)• 


LL 




dK(LL - K) 
(LL - K)^ + 


F (K,J - LL + K - N) 


and 


(4-9) 


E(J) = 


LL“ 


2N^k^ 


■/i ^ 


f(ll + 1, Jj) 



dK K 


+ 2N^k^ 
k(k^ + 4N^)^ 


- 1 


F (K, J - K - 2N , (4-10) 


with 


J, = J - LL - N for large and moderate values of © 
small values of 


© 


and J, = J - LL for 


For a numerical computation scheme, I, X, and J are taken as integers; 
while LL is always maintaihed as an even integer. The integrals in K can 
be represented ^s a sum with the integrand evaluated at odd values of 
K = 1, 3, 5 ... (LL -1), so that dK = 2. However, this procedure had to 
be modified to account for the rapid variation in the kernels of some of 
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the integrals. Thus, in Equation (4-8), the conventional way of interpreting 
the second integral would be as 


LL-1 


^ ^ 1 - |k - l| [(K - 1)2 + 6N2k2jy'[(K - 1)2 + 4N2k2)] | 


N=1,3 


X f(K, J - |k - l| - 2N) 
However , the kernel 1 1 - | K - 1 1 j^(K - 


1)2 + 6N2k2j^jjK - 


1)2 + 4N2k2j 


is a steeply varying function of K with a symmetry around K = I, where 
it has a maximum value of 1. The above summation is hence an overestima- 
tion. A more realistic interpretation would be to evaluate the kernal at 
]k - 1 1 = i rather than K = I, and hence represent the integral as equal 


to 


LL-1 


JL V 

NR 2 ^ 


1 - K - I 


[<K- 


1)2 + 6N2k2j j ^ 


(K - 1)2 + 4N2k2] 


K=l,3. 

K4I 


X F(K, J - Ik - I - 2N) + — 
' Nk 


1 - 


(l_j_24Nf)^) 

(l + 16 N 2 k 2 ) ^ 


F(I, J - 2N) 


Further, in Equation (4-10) , the terms 


k 2 + 2N2k2 


(k2 + 4N2k2] 


-K 


and F(K, J - K - 2N) , representing the integrand of the second term, are 
both functions that rapidly drop to negligible values as K increases. 
Designating the kernel 


N2k2 


r 2 + 2N2k2 

fK2 + 4N2k2')^ 


- K 
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as P(K), and F(K, J - K - 2N) as g(K) (only space variation is of concern), 
the integral under examination, viz. 


LL 

Jds P (s) g(s) , 


is interpreted as, 
LL-1 K+1 

I I 

K— 1 , 3 , . , . (K— 1) 


ds P(s) [g(K) + g-(K) + g’» (K) + 


-] 


where g(s) is expanded as a Taylor's series around S = K and g' (K) , etc., 
are the derivatives with respect to S evaluated at K. Ignoring g" (K) " and 
higher derivatives, and considering g* (K) only at K = 1, the integral 
under consideration is representable as 



F(K, J - K - 2N) + |f(3, J - 3 - 2N) - F(I, J - 1 - 2N)J 





X 
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'^hus. Equations (4-8) through (4-10) take the following form; 
tan I 


F(I, J) = - 


f 


J 


d0 sin^9 



LL-1 


^ Nk 


K=l,3,— 

K^I 




|k - l| [(K - I) ^ + 6N^k^] 


(K - I)^ + 4N^k^]^ 


F{K, a - !«■ - ll - 2N) 


Nk 


1 - 


(1 + 24N^k^) 


(l + ISN^k^)^-* 


F(I, 


(LL - 1) 
2Nk 


(LL - + 2N^)c^ 


(LL - I) (LL - I)^ + 4N^=k 


2 , /„>,2,.2j% 


- 1 


X F(LL +1, J-LL+I-N) 


F(LL + 1, J) = e VJ 


1 - 


e 


1 + 


N^k^ 


LL 


LL-1 


+ 4N^k^ ^ 


(LL - 


k=i,3,-[(ll - k) 


F (K, J - LL + K - N) 


tan^e \^ 
4N^k2 ) 


J - 2N) 


(4-11) 


K) 

2 + N=^k^]2 


(4-12) 
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and 


E(J) 


•KS) (■ • (‘ • ) 


f(ll + 1, Jj) 


+ - {f(3, J - 3 - 2N) -F(1, J - 1 - 2N) 





K=l,3,— 


(K - 1) / (K - 1) 


2Nk V ^2,^2 




+ 4 



F(K, J - K - 2N) 


(4-13) 


where Jj = J - LL - N for large and itraderate values while Jj = J - LL 

for small values of 


(i)- 


(i)' 


Equations (4-11) through (4-13) can be easily solved numerically by 
the process of "stepping forward in time." The first integral in Equation 
(4-11) is evaluated using a five-point Gaussian quadrature scheme. Due to 
the numerical finite difference quadratures [F(I, J) , F (LL + 1, J) , and 
E(j)] in the limit J -> <» do not tend to 1 as they should. They are hence 
forced to tend to 1 by using correction factors: HF(I), the correction 

factor for F(I, J) , multiplies all terms on the right side of Equation (4-11) 
except the first integral: HB and HE, the correction factors for F (LL +1, J) 
and E(J), respectively, multiply Equations (4-12) and (4-13), respectively, 
with 


_ -1 
tan 

i 


2Nk 


d9 sin^0(l 


L _ tan"e \^ 
V ' 4N2k^/ 


HF(I) 


denominator 
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where 


denominator ; 


(LL - I) 

(lL - i)^ + 2N^k^ 


. j_ 

1 — 

(l + 24N^k^) 

2Nk 

(LL - I) ^LL - I)^ + 4N^k^ 

I* J 

Nk 

1 "• 

(l + 16N^kO^- 


LL-1 



K=l,3,— 

K=I 


(■ 


lil 


- l| [k - I)^ + 6N^k^] 
[k - I) ^ + 4N^k^] ^ 


(4-14) 


HB = 




LL^ + N^k^ 


+ 4N^k^ 


-1 


(LL - K) 


K=l,3,— 


[(LL - K)^ + N^k^J^ 


(4-15) 


and 


HE = 

/ 

denominator 

where 

denominator : 


1 /LL\^ (-/ 2N^k2x r 4N^k^-l"l V 


4N^k^ ' 


LL-1 


K=l,3,— 


( k + l) 

2Nk 




% (k - l) 


,2 


2Nk 


r-si") • •] 


(4-16) 
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The time constant is determined as that value of J at which E(J) 
exceeds (1 - 1/e) , Two computer programs were formulated; one to solve 
Equations (4-11) through (4-13) for large and moderate values of LL/Nk 
ranging from 24 to 1, while the second was used to solve the same 
equations for small values of LL/Nk ranging from 1 to 1/64. It is 
worthy of note that using finer and finer meshes by increasing N was 

not found to change results at all (increasing N is equivalent to taking 

I 

finer meshes in space as also finer increments in time) . 

Figure 1-2 shows a plot of the time constant (upper bound) vs LL/Nk. 

Figure 1-3 shows an amplified plot of Figure 1-2 for small values of 

(LL/Nk) . 

Lower Bound to the Time Constant 

The relevant equations are again (4-4) through (4-6) , and the first 
two assximptions made for the upper bound are made here again. Assumption 3 
for the upper bound case is still used, but the "averages" are deliberately 
underestimated. This is equivalent to \anderestimating the time of flights. 
Thus, in Equation (4-4) we represent the average over 6 of 

- z)^ + 2a^(l + cos6^ 
by I z ' - z I cind 

- z^^ + a^ + r^ + 2ar cossj 

at r = 0 (due to assumption 2) by (Z. - z) . In Equation (4-5) we represent 

- z') ^ + a^ + r^ + 2ar cosoj 

at r = 0 (by assumption 2) by (£ - z'). Finally, in Equation (4-6), we 
represent 

j^z^ + a^ + r^ + 2ar cossj 

and 

+ r^ + r'^ + 2rr' coss]^ 

at r' =0 (by assumption 2) averaged over 0 and r by z and 1, respectively. 
Scaling is done as in Equation (4-7) with k now set equal to 1. In finite 
difference form, the resultant equations (after amendments of some of the 
integrand as in the case of the upper bound) are: 
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Fig. 1-3. Variation of time constant with small length to radius ratio. 
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F(I, J) 


F(LL + 


E(J) = 


+ 


.+ 


— 1 

tan ~ 




J,. Il , tane ..\* , X y 


LL-1 


_ |k - l| (k - l) 


^ + 6N^ 


2T2 


K+I 


|^(k - l)^ + 4N^J 


F(K, J - |k - l() 


N 


1 _ (l + 24N^) 
(l + 16 N^) 


F(I, J) 


(LL - I) 
2N 


(LL - I)^ + 2N^ 


(LL-I) |^{LL - I)^ + 4N^J' 


- F(LL + 1, J - LL + I) 


(4-17) 


1, J) = 


e VJ 



LL-1 


+ 4N 


s 


(LL - K) 


K^l.X,- [<XX.- 


F(K, J - LL + K) 


(4-18) 


-(— V 

2 ' N / 


I 2N^ /, 

1 + - (1 + 

LL^ V 


4N^\ ^ 

LL^/ 


F(LL + 1, J - LL) 


3 


F(3, J - 3) - F(1 , J - 1) 




^ 1 
+ 2N 


LL-1 

I 

K=l,3,- 
\ 

2K , , 

— }F(K, J - K) 


(K + 1) 
2N 


(K 


. 4 ]^ 


(K - 1) 
2N 


(K - 1) 


+ 4 


(4-19) 
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Correction factors HF(I), HE, HB from Equations (4-14) through 
(4-16) are used again. Another computer program was evolved to solve 
Equations (4-17) through (4-19) , with the time constant being determined 
as that value of J at which E(J) exceeds (1 - 1/e). Figure 1-2 shows 
the plot of the time constant vs LL/Nk. 

It bears observation that the time constant obtained is truly the 
lower bound since the length of flight, and hence the time of flight 
used, is the minimum possible (within the framework of assumptions 1 and 
2). In fact, the time of flight between two elements ^at a speed 
is assumed to be equal to that of a particle moving parallel to the z axis. 
We hence have the anamoly that the speed of particles in the' radial 
direction is infinite. However, we have achieved bur objective — that 
of finding an upper cuid a lower boxind to the time constant, although as 
evidenced from Figure 4-2, the span between the two bounds is rather large. 

Conclusions 

The plot of the time constant vs 2./a shows that the time constant tends 
to a constant as i/a increases. This is directly attributable to the 
assiamption of diffuse reflection from the wall wherein each wall element 
reflects to the exit a fraction of the flux coming to the element/ whose 
value depends on the solid angle made by the wall element with the exit. 
Thus, while the time constant increases with i/a for small i/a, after a 
certain stage adding new wall elements (i.e. , increasing the cylinder 
length) is not going to increase the time constant, since not only does 
the solid angle subtended by the wall element to the exit decrease rapidly 
(as axial distance from exit increases), but also the flux to the element 
(and hence flux from the element) is small in the time scales considered. 

It is obvious that, if, the walls were specularly reflecting, the time 
constant will progressively increase with i/a. Most svirfaces are partly 
diffuse and partly specular reflecting, being more diffuse than specular 
in nature. It is hence possible to anticipate results similar to our 
results for real surfaces; however, the exact value to which the time 
constant tends, depends very much on the nature of the cylinder surface. 
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APPLICATION OF SUDDEN FREEZE MODEL TO POROUS MEDIA 


Formulation and Assumptions 


The degree of sophistication necessary to analyzing the flow through 
porous media with vacuum as one boundary condition is evaluated. It is 
recognised that the f low , . in general , would contain an initial region of 
continuvun flow that would undergo a gradual transition to free molecular 
flow. However, if the initial region of continuum flow prevails over a 
sufficiently large distance from the probe, the continuum equations can 
be used with very little loss of accuracy. In an effort to obtain rough 
estimates of the length of this continuum flow, a sudden freeze model is 
proposed. Such a flow model is applied to various models of the pOrous 
media. 

The basic assumptions made here are : 

1. The temperature T of the flow is constant. 

2. The flow which is initially continuum undergoes transition to free 
molecular flow abruptly at a section termed "the freeze section." 
Properties at this section are described by the subscript "f." 

The freeze section is identified as one wherein the mean free path, 
X, is equal to the pore radius, a. 

First a one-dimensional flow through a slab of finite thickness fi, is 
considered. The medium is visualised as being made up of a bundle of 
straight, capillary tubes of constant radius, a . Let the pressure at the 
entrance of the tube (i.e., at the probe) be p^ and the pressure at the other 
end of tube be equal to 0. Let the freeze section at which the mean free 
path X = a, have a pressure p^ and be at a length from entrance. 


In the continuum section, the flow is a poiseuille flow, 
mass flux through each tube per unit time is 

4 


Hence the 


Q = - 9 


TTa 

8y 


dx 


(4-20) 


where y is the viscosity, dp/dx the pressure gradient along the tube axis 
cUid X the coordinate along the tube axis. Integrating (4-20) with respect 



to X from the entrance to the freeze section, we have (since p 
Q is a constant) 




Assuming Pf have 


O = 


2 4 

p TTa 
1 

16y RT 


In the free molecule region 

^ ^ ‘SX 2 

Q = - C — a TTa 


where C is a constant of 0(1) and V 


_ / 8RT\^ 


\ IT / 


Integrating (4-23) between the freeze section and the tube 


have 


Q = 


C TTa^ ^f 


V 


( «- %) 


Since at freeze section A = = 


2y 


f /P 


(^) 

\RT / 


= a. 


hence substituting for p^ in (4-24) and equating (4-22) and (4- 
we have 


SL 


1 + 8C 


(^)T' 


where 


A. = 


2y 


H-) 

\RT / 


= mean free path at tube entrance (i.e., probe). 


= 9RT and 


(4-21) 


(4-22) 


(4-23) 


exit, we 


(4-24) 


24) for Q, 


(4-25) 
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Next, a one-dimension spherically symmetric flow through a semi-infinite 

\ 

medium' is considered. Let r^ be the probe radius, and p^ and X, the 
pressure and mean free path at the probe. (Strictly^ the end of the probe 
will probably not be a hemisphere, so r^ should merely be regarded as a 
typical length seal for the tip of the probe) . Let r^ be the freeze 
radius and p^ the pressure at the freeze radius . Let the pressure at 
infinity be zero. The porous medium is visualised as an assembly of 
isotropic, randomly oriented, straight cylindrical pores of constant 
radius "a", connected to one another at the ends. Several pores may start 
or finish at these end points. 

In the continuum region the mass flux Q is given by Darcy's Law as 

Q = - 9 !^^ 2TT r^ (4-26) 

U dr 


where K is the permeability and r is the radius. , 

Integrating (4-26) between r and we have (noting 9 = ^ and that Q 
is a constant) 



Q = 


TTK 

yRT 




(4-27) 


We now consider the free molecular region. Consider a point at radius r. 
The pores being randomly oriented, the number of tubes at radius r with 
angular orientations (spherical coordinates r, 0 and (j)) between 0 and 0 + d0 
and 0 and 0 + d0 is 



sin9 d9 d0 
4 tt 



sin9 d9 d0. 


where a is the porosity of the porous medium. 
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Since, in the free molecule limit, the mass flux through each tube is 



COS0 §2. ^ 

dr 


hence the mass flux over a radius r in the outward direction (along 


increasing r) is 




CTTa^ 

V 


COS0 ll) 



dr 


(4-28) 


Integrating (4-28) (and noting that Q is a constant) between r^ and infinity, 
we have 


ar. 


Q - 7T c a 


(4-29) 


Equating (4-27) and (4-29) , assuming p » and using the fact that 


^f = a = 


2yRT 

VP, 


we have 


— = 1 + 


2K 


c a Xj 



where we assume that K = M a^ with M, a constant of 0(1). 


(4-30) 
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Finally a one -dimensional flow through a slab of thickness £ is 
considered, using the above random pore orientation model. The notations 
used are the same as for the case of parallel capillary tube with £ ^ now 
denoting the section at which flow freezes. 


In the continuum region, the mass flux Q per unit are given by Darcy's 
law as 



9K ^ 
y dx 



where x is the coordinate along the line joining the end faces of the slab. 

Integrating the above equation between x = 0 (from probe) to the freeze 
section x = S,^ aind using 9 = ^ we have. 


K (p;-<) 

2yRT £ 


K 

2yRT 



on assuming Pj ^ ^ P 3 • 


(4-31) 


Let us now consider the free molecule flow. Here at a section, x, 
the total number of tubes with angular orientations between 0 and 
9 + d9 , <() and 4 > + d(|) , per unit area is 



sin9 

4tt 


d9 dcf) . 


Mass flux through each tube 


= - C 


dx TTa" 


cos9 


Hence total mass flux Q per unit area at x is 


J 27T 


9=0 <|)=0 


a c a dp 
^ ^ ■ 


(4-32) 
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Integrating (4-32) with respect to x from x = 2,^ to x = jj, (and noting Q 
is a constant) we have 


g C a 

4V “ ^f) 


Equating (4-31) 


and (4-33) , 


and noting that 


2yRT 

Va 


we have 


(4-33) 



(4-34) 


where 


A = 
1 


2M 


v(M 

V RT / 


= mean free path at tube entrance (i.e., probe) 
while we assume K = Ma^ where M is a constant of 0(1) . 

CONCLUSIONS 


To s\immarize the results obtained, we have: For the one-dimensional 

slab flow using both the porous medium models, we find that the ratio of 
the freeze length to the total slab thickness is ^1 + 0(1) ^X^/a^ J \ 

where X, is the mean free path of the gas at the probe entrance and a the 
pore radius. In the spherically semi-infinite flow, the ratio of the 
freeze radius to the probe radius is estimated as + OU) 

— 2 —3 

On the basis of a moon grain size of 10 to 10 cm, the pore radius 

-3 -** 

can be roughly gauged to be 10 to 10 cm. If, for example, nitrogen 
is pumped from the probe to the lunar surface at a pressure of about 1 
atmosphere and normal (15° C) temperature, then X^ » a, and transition 
will occur far from probe. Hence the continuum equations can be safely 
used. However for rocks of smaller pore size (£ 10 ^ cms) X^ and a are 
of the same magnitude and transition will occur near the probe — a detailed 
analysis of the transition flow is called for in such a case. 
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APPENDIX 

1. Formula for entrance to cylinder flux C (z, t) 

c 

The initial value for the distribution function f at time t = 0 is 

f = exp I 1 6 (z) , where ^ -is the molecular velocity 

and 6 (z) is a function such that 

6(z) = 1 for z < 0 

= 0 for z > 0 . 

Solving the unsteady Boltzmann equation in the free molecular limit for 
the above initial value, we have 

/-§ § 

S' S' « - 7 ^ 

(2tt RtJ2 



If n be the unit normal at the wall at z, the number flux from 
the entrance to an annulus of unit area located at the wall at z, at 
time t is 



Using a spherical coordinate system fixed to the wall as shown 
in Figure 1-Al, we have 

? • n = 5 sin0 cos4> . 
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Employing the ray-tracing technique gives 



The geometry for the entrance to cylinder flux is shown in Figure 1-Al. 



Fig. 1-Al. Geometry for entrance to cylinder flux. 
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2. Formula for cylinder to cylinder flux ^ (z, t) 

cc 

By virtue of the basic assumptions and definition of t) , the 

distribution function of the molecules emitted from the cylindrical wall 
at time t and axial position z is given by 


f 

c 


(2, t) 



The geometry for the cylinder to cylinder flux is shown in Figure 1-A2. 


Consider two elements of area dA and dA' at axial positions z and z' , 
respectively, separated by a distance Sj,- Let <)) and (})' be the angle made 
by the line joining the two area elements with the normals at dA and dA', 
respectively. Then 

= (z' - z)2 + 2a^ (1 + COS0) 

2 

, . , a(l + COS0) 

coscp = coscp' = : g 

d A' = a d0 dz* . 


S 

Molecules in the velocity range ^ and ^ + d^ leaving dA' at _ t - 
can arrive at dA at time t. Hence the number of molecules with velocity 
in range ^ cind ^ + d^ leaving dA* that arrive at dA at time t 


= dA d^ (z, t) 
cc 






/ t 



COS(l) cos(j)' 



dA dA' 




Fig. 1-A2. Geometry for cylinder to cylinder flux • 
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or the total number of molecules in all speed ranges that arrive at time t 
on an unit area annulus at z per unit time from all parts of the cylindrical 
surface is 


■// 


~ I j cos(|) cos4>’ I d^ 5"* f^ 

A- 


^c(^ ' " F^\ 


= ^ J d0(l + cos8)^y^ 


dz ' 


J [(z* - 
0 *- 


z)2 + 2a^(l + cos0)J 


00 

/ 


X I dc c^ e ^ ^ iz’ t - 

• ’c * 


f 


[(Z* - 


z)2 + 2a^(l + COS0)1^ 




3. Formula for back to cylinder flux C. (z, t) 

DC 


Pursuing the same line as in the C^^(z, t) derivation and referring to 
Figure 1-A3 , we have here 


§, t) = §» 




Sg s= (A - z)^ + a^ + r^ + 2ar cos0 
cos(J) = (a + r cos0)/S^ 

cos<p’ = (Z - z)/Sj 

and 


dA' = rdr' d0 . 






1-32 


Hence, 




, t) =J^ cos^ costp J ^b(^' §' * ■ *^) 


IT “ 00 ^ 

— fdQ fdr r f ^ ^ (g, - z) (a + r cos6) 

^ J J J f(£ - z) ^ + 2ar cosSl^ 

0 0 0 L J 

[(g, - z)^ > a^ + + 2a'r cossj 




(2%)- 


To get the upper and lower bound (of the time constant) approximation, 
using assun^tion 1 . and 2 . , we have 


r, t - 


[ 


(i - z) ^ -f + 2ar cos^ 


J . U - 2)^ * a"] 


1 


SO that 


TT a 





dr r(g, - z) (a + r cos9) 

|^( Jt - z) ^ + a^ + r^ +2ar cosbJ ^ 


/ 


, 3 -“C 

dc c e 


= -( A - 

2a 


(g, - z)^ + 2a' 


(£ - z) [(i - z)^ + 4a^ 


- 1 


X t 


[(& - z)^ + aj ^] 
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4 . Formula for cylinder to back flux ^ , (r, t) 

cb 

Referring to Figure 1-A4, we have here 

= (£ - z ’ ) ^ + a^ + r^ + 2ar cos0 
cos(j)' = (il - z' )/S 

4 . 

cos<t> = (a + r cos9)/S 
dA' = a d9 dz* 



To get the upper and lower bound (of the time constant) approximation, 
using assunption 1. and 2., we have 



0 
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5. Formula for entrance to back flux ? . (r, t) 

eb 

The number flux from entrance to the back end at time t at radius r 
per unit area is. 


'eb*'- 


ff f ^“§(5 • 2)fe 


where n is the unit normal to dA . 

Proceeding on lines similar to the calculation of C (z, t) and 

ec . 

referring to Figure 1-A5 we have 


a+r r^-a^+jl^tan^0 
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To get the upper and lower bound approximations (for the time 
constant), using approximations 1,. and 2., we have 





0 0 J, 


tcos6 



Formula for cylinder to exit flux e^(t) 


Referring to Figure 1-A6, we have here. 


C (z, t) 2 (tt) 

f(r, Cf t) = f (z, t) = 

- - c - 2 

/ Nfc / \ 2 


(2ET„)(2n RT„)‘ 


= z^ + a^ + r^ + 2ar cos0 

6 

cos(() = (a + r cos0)/S 

6 

dA = 2TTa dz and dA' = rdr d0 . 


roico 



y/imimmimiiii// 


1 

1 

a 

Y dA 

i, 



W7mm 

-It 



Fig. 1-A6. Geometry for cylinder to 



it flux . 





since 


, dA dA- . COS.I. cosj,- J / E, t - 

' s^ '\ ■ s/ 


or 


e (t) 
c 


‘f^'f 

A' A 


dA cos(|) cos4> 


-/ dEE= f^(z. 


= 8a 


JL a TT 


ao^a + r coso 




+ a^ + + 2ar cosej^ 


•/ 


- -c 
X f dc c e 5 


t (z^ + a^ + + 2ar coss)^ ! 

" ■ J 


The geometry for the cylinder to exit flux is shown in Figure 1-A6 


7 . Formula for back to exit flux e, (t) 

b 


Referring to Figure 1-A7, we have here. 


f(r, t) = fj^(r, §, t) = 


Cj^(r, t) 2 (tt) 




\ 


3 exp 

2 \ 2RT 


w 


+ r^ + r'^ + 2rr' cos0 

6 


cos(|> = cos<t>' = S-/S 


dA = 2TTr dr 


and dA' = r' dr' d0 
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Since 


de^ = dA dA- dC 


hence integrating we obtain 

tr a 


e^(t) 


= 8 r dr r r de f ^ 

J J J + r^ + r’^ + 2rr' coser 
0 0 0 > J 


00 

■/ 


2 > . Vi + r + r'" + 2rr' cos 

3:/t - — ; 

dc c e b „ 


se ) 


:(2RT ) 
' w' 


To get the lower eUid upper bound approximations (for the time constant), 
using assunptions 1 and 2, we have 


r, t 


CM 

+ 

CM 

+ 

CM 

c»? 

■ 2rr* cosb) 


L . ' 


C(2RT„) 

1* J 

^ f w - 

L 1 

H»)* 



so that 


e (t) = 4i 
b 


IT a a 

f' ^ dr r 0, t 

^ / d0 / dr r' / 

J J J [i^ + r^ + r'^ + 

0 0.0 


dr r 

Fo. 

[i^ + r^ + r 

+ 2rr' cosg]^ 


The geometry for the back to exit flux is shown in Figure 1-A7. 
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SYMBOLS 


a 

c 

e (t) 

e. (t) 
1 

E(J) 


F. 

1 

F(I. J) 


F(LL + 1, J) 

HB 

HE 

HF(I) 

I 

J 

k 

K 

s, 

h 

m 

M 


N 


n. 

1 


Pi 


Q 


radius of pore — dead end, open. 

scaled modecular speed; constant of 0(1) 

total outflux from dead end pore at time t 

outflux (from dead end pore) from i at time t 

scaled outflux from dead end pore 

distribution function of molecules coming from i 

scaled cylindrical wall flux 

scaled flux from the back end of the dead end pore. 

correction factor for the back end flux 

correction factor for the outflux 

correction factor for the cylindrical wall flux 

scaled distance along dead end pore axis 

scaled time (from opening of the cylinder entrance) 

constcint 

permeability 

length of the dead end pore; length of porous medium slab 
length of freeze section from, probe 

square root of the ratio of wall temperature to temperature 
of gas outside the dead end pore 

constant of 0(1) — relates permeability and the square of 
the pore radius 

an integer 

number density of molecules coming from i. 
pressure of gas at i 
mass flux 



r 

r. 

i 

R 

t 

T 

T. 

1 

V 


z 

a 



X. 

X 

y 

K 

9 


radial distance from dead end pore axis; radial distance 
of spherically symmetric flow. 

radial distance of i in spherically symmetric flow 
gas constant 

time from Opening of dead end pore entrance 

scale for time, temperature 

teir^jerature of molecules coming from i. 

most probable molecular velocity 

axial distance along dead end pore axis 

(alpha) porosity 

(zeta) molecular flux from i 

(zeta) modecular flux from i to j 

(lambda) means free path at i 

(mu) viscosity 

(xi) molecular velocity 

(rho) density 


SUBSCRIPT 

b 

c 

e 

f 

w 

i 


back end 
cylindrical wall 
entrcince 
freeze section 
wall 

probe entrance 



Chapter 2. EXPERIMENTAL STUDIES 
F. C. Hurlbut, C. R. Jih, and P. A. Witherspoon 

INTRODUCTION 

The underlying rationale for undertaking studies of fluid flows 
in porous media under rarefied gas flow conditions has been to supply 
the en^irical basis for theory necessary to the design and understanding 
of a permeability probe device for in situ experimentation at the lunar 
surface. A second objective has been to provide the practical experience 
in such experimentation necessary to permit a sound, efficient, and work- 
able design of such a probe. 

The outlines of our attack, both theoretical and experimental, have 
been summarized in the 1969 Final Report, Vol. IV of IV, "Studies on 
Conductivity of Lunar Surface Materials," by Katz, Willis, and Witherspoon 
and remain very little changed to this date. In that reiport our current 
state of knowledge was described, preliminary concepts of probe design 
were discussed, and directions of analysis were indicated. 

The present report is confined to the description of the ongoing 
experimental program and to a presentation and discussion of preliminary 
observations. It should be understood as a record of work in progress 
and is to be taken, together with the 1970 Final Report, "Studies of 
Fluid Conductivity of Lunar Surface Materials —Theoretical Studies," 
by Raghuraman and Willis, as a representation of our progress in the 
fiscal year 1969-1970. 

BACKGROUND 

The flow of gases through porous media has received a moderate 
amount of attention over the years as, for example, the flows of low- 
density gases connected with the problems of catalytic beds or with 
those of transport phenomena at permeable barriers. The words "low 
density" refer here to conditions under which the Knudsen number, based 
on pore size, is greater than 't 1/100. One may note that such low- 



density flows might well occur at pressures above or below 1 atmosphere 
for rocks within the ordinary range of pore size. Above the limit of 
low density cited, the flow of gases in porous media may be treated by 
the empirical continuxim methods which have been found to be successful. 
Prior studies of low-density flows have been confined to certain semi- 
continuum models or to the assumption that the flows are entirely free 
molecule in character. Such models imply that the density gradients 
are everywhere small, a constraint which cannot be applied in general 
to flows in porous media whose natural environment is, and has been for 
a very long time , a vacuum . , 

Related studies of the flow through capillaries have been more 
widely conducted, and it would be in connection with these somewhat 
simpler flows that one would hope to see the development of theoretical 
models for the transition from continuum flows to the free molecule 
regime. Such models would provide a valuable base for modeling the 
porous medium. However, we again find that nearly all theoretical work 
has confined itself to conditions where the density gradients are small 
so that the gas remains within a particular regime of flows throughout 
the capillary. Work relating to larger density gradients has been con- 
ducted by interpolation and fitting but without a rigorous basis in the 
kinetic theory. Experimental work on capillary flows has been conducted 
under conditions appropriate to the theory with few exceptions and in 
these latter cases no examination has been made of the details of the 
transition from continuum flow to free molecule flow. 

With these limitations of available information in evidence it was 
determined to undertake direct measurements of porous medium permeability 
under low-density conditions as the most efficient route to the design 
and understanding of an in situ permeability probe for lunar materials. 

It was determined that initial investigations should be of one-dimensional 
flows through homogeneous, simulated rock samples having a range of 
permeabilities. Use would be made of the pumping system- associated with 
the existing rarefied gas wind tunnel , and it was also planned that 
advantage would be taken of the technology and experience of the Rarefied 
Gas Laboratory.* The program of design, construction, and measurement 
has proceeded well, but not as rapidly as planned, so that to this time 

*U. C. Division of Aeronautical Sciences (Mechanical Engineering) 
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only the first phases of the measurement program have been completed, 
the next sections details of the permeability apparatus are given. 


DESIGN AND CONSTRUCTION OF APPARATUS 
Introduction 

As in the proposed permeability probe , gas from a source at moderate 
pressures flows into the porous specimen toward a sink at low pressures. 

If the Knudsen number of the flow is initially of order 1 or smaller, 
the flow will inevitably transform to a free molecule flow within the 
specimen. The character of the transition, as determined by the measured 
pressures at various distances from the source, will permit a calculation 
of permeability and possibly pore size and configuration when suitaible 
theory becomes available, The experimental apparatus required for the 
investigation of one-dimensional flows within the above conceptual 
framework consists of a gas source and flow metering system, a specimen 
chamber with pressure taps distributed along its length, a pressure 
transducer and metering system, a high capacity vacuum pump, and the 
necessary valves, ancillary gauges, and pumps. A detailed description 
follows . 

Description of the Apparatus 

It was a basic objective of the design that it should permit the 
detailed examination of pressures as a function of position along a one 
dimensional flow through a porous specimen. It was determined that the 
specimen should consist of up to 10 cylindrical slabs, each of thickness 
to 1 inch and diameter 2.5 inches, permitting pressure measurements to 
be made at discrete intervals by sampling the space between slabs. The 
arrangement is shown in Figure 2-1, a dimensioned assembly drawing of 
the equipment. The specimen chambers, shown in greater detail (Figure 
2-2) , consist of 2 flanged cylinders of stainless steel each with pro- 
visions for 5 segments of specimen. Each specimen segment consists of a 
plexiglass ring within which is cast the porous material . A seal between 
the plexiglass ring and the inner wall of the specimen chamber is arranged 
by an "0" ring set in a groove in the chamber wall. Pressure taps with 
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Fig. 2-1. System assembly. 







pressure leads of 1/4" diameter stainless tubing welded in place are 
provided between each specimen position. Spacer rings between each 
specimen maintain the correct position. 

Uniform entry conditions over each slab face are established by 
virtue of the high flow conductance of the large gap between slabs as 
compared with the lower conductance of the slab material. Thus a 
segmented ideal one -dimensional flow is permitted. Note that the 
number of slabs may be varied from 1 to 10 and that the thickness of 
each slab may be arbitrarily determined up to 1 inch. 

Each pressure tap is connected via a valve to a central manifold 
and that manifold is connected through 3/4-inch copper tubing and a 
quarter swing valve to a pressure transducer. Ample conductance is 
provided to permit degassing the specimens and to make possible a 
sufficiently short gauge response time. The pressure transducer is 
an MKS diaphragm gauge having a maximum differential pressure range of 
30 Torr. This device was selected for its well-known accuracy, stability, 
and insensitivity to gas composition. Since it is a differential pressure 
gauge it must be connected to a reference vacuum system, details of which 
are shown in Figures 2-1 and 2-3. Note that a bypass valve is provided 
for establishing the zero of the gauge and to permit evacuation of the 
manifold. 

At the downstream end of the specimen chambers is a 6-inch vacuum 
gate valve and beyond that, the main manifold of the rarefied gas wind 
tunnel. ' The pumps associated with the wind tunnel flow system have the 
capacity to maintain the downstream end of the permeability apparatus 
at 1 to 2 microns Hg for any realistic flow within the porous samples. 

The gas supply system is also shown in Figure 2-1. This system 
consists of high and intermediate pressure regulators, appropriate shut- 
off valves, and a gas-service regulator followed by a system of 5 viscosity 

- 3 

type flow raters covering a range of flow rates to 3 x 10 cc/sec. 

For lower rates of gas flow the film-capillary, positive displacement 
method will be used. A needle valve between the flow metering system 
and the first specimen chamber serves to regulate the flow rate. 
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Fig. 2-3. Pressure metering system. 
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An overall view of the permeability apparatus is shown in the 
photograph, Figure 2-4. 

SPECIMEN PREPARATION 

A number of options exist for the preparation of porous samples. 
Among these are the aggregation of sand particles using wax or resin 
binders; the casting of concretes and artificial stones; the sintering 
of metallic beads, chips, filings, or fibers; the sintering of ceramic 
materials in fibre, bead, or rod form; and the cutting of natural rocks, 
particularly those of volcanic origin. To be suitable for low pressure 
measurements, the porous media must be free of organic materials having 
vapor pressures in the micron range. In order to insure an extended 
region of transition flow under pressure conditions appropriate to 
this experiment, the samples should have much greater permeability than 
ordinarily found in natural rock. 

The various constraints of the present program favored the con- 
struction of sintered materials, preferably ceramics. However, for 
initial performance testing it was decided that cast concrete specimens 
would serve and that these could be constructed using materials and 
technology readily available in the Civil Engineering Laboratories. 
Accordingly, three sets of cast concrete samples were prepared, careful 
attention being paid to mixing and uniformity of casting procedure. The 
conposition of these samples is shown in Table 2-1. 

Table 2-1 


Saitple 

Sand 

(gm) 

Cement 
^ (gm) 

Water 

(gm) 

1 

2000 

100 

100 

2 

2000 

500 

250 

3 

2000 

80Q 

400 


In each set the concrete was cast into 12 to 16 plexiglass ring 
forms of 1-inch depth. Upon curing and drying the samples, the perme- 
ability for air was measured at atmospheric pressure and above using 
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Fig. 2-4. View of the permeability apparatus. 




a standard permeability apparatus. Set 3 was rejected immediately 
as too impermeable, and the 10 slabs most uniform in permeability were 
selected from each of sets 1 cind 2. The results of these tests are 
shown in Table 2-2. 

Table 2-2 


Slab 

No. 

Permeability 

K(cm^) 

(X 10~*^) 

Slab 

No. 

Permeability 

K(cm^) 

(X lo"'^) 

1-1 

9.3 

2-2 

4.5 

1-2 

9.3 

2-3 

6.8 

.1-3 

9.3 

2-5 

' 4.5 

1-4 

9.3 

2-7 

4.8 

1-5 

9.1 

2-8 

7.5 

1-6 

9.1 

2-9 

4.8 

1-7 

9.1 

2-10 

5.2 

1-8 

9.1 

2-11 

5.8 

1-9 

9.5 

2-13 

5.2 

1-10 

9.5 

2-16 

5.4 

It may be 

noted that Set No. 

1 is both 

more uniform and 


permeable than Set No. 2. The permeabilities are within the range 
of the more porous natural rocks. 


PRELIMINARY OBSERVATIONS 

Preliminary measurements were made on a set of 4 slabs (No.'s 
2-2, 2-3, 2-10 and 2-11) to gain operational experience with the 
instruments and to develop a physical sense for the ■ appropriate 
permeability ranges of the next generation of specimens. The gas was 
nitrogen. The tunnel downstrecim pressure and the M.K.S. gauge reference 



pressure were at 1 micron or below. Operational experience was 
obtained , although no useful quantitative information has resulted to 
this point, owing to the low permeability of the present specimens. 
Certain conclusions may be drawn which are summarized as follows: 

1. It must be recalled that the objective of the system 
design is to permit the study of the regime of transition 
flows in porous media. This is accomplished by extending 
the region of transition in physical space over slabs 
which are somewhat less in thickness than that region. It 
is implied that specimens of very high permeability must 
be eitployed. 

2. The flow conductance of the specimens must be sufficiently 
great that the pressure taps open into, effectively, an 
unlimited reservoir of gas at the measured pressure. 

t . 

Thinner samples and greater permeabilities will improve 
conditions in this regard. 

3. In all regards the apparatus behaved well and appears to 
have the capability of giving results of the quality desired. 

CONTINUING PROGRAM 

Within the next few months it will be our objective to complete 
measurements enabling the description of transition flows in porous 
media. Interpretation of these results will require the independent 
characterization of the medium in terms of pore size and configuration. 
Such characterization will be accomplished by a combination of optical 
and displacement methods and by a knowledge of the size and configura- 
tion of particles (beads, rods, etc.) used for the preparation of each 
sample. Materials of various descriptions will be formed into specimen 
slabs, with sintering being viewed as the most promising technique at 
this time. Thus, one may summarize by stating that the next phase of 
our activity will consist of four essential parts: 1) the preparation 

of suitable samples, 2) the physical characterization of these sanples, 
3) the measurement of flow characteristics, and 4) the interpretation 
of results. 



